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nuclear temperature constant, MeV -1 
energy level density parameter, MeV -1 
energy, MeV 
electromagnetic 

giant dipole resonance central energy, MeV 

photonuclear reaction threshold for a particular species X, MeV 

photon energy, MeV 

electric multipole of multipolarity l 

fraction of direct, nonstatist ical emission 

one-particle level density at Fermi energy, MeV -1 

energy independent branching ratio for species X due to 7 rZ nuclear 
excitation 

energy dependent branching ratio for species X due to nl nuclear 
excitation 

Boltzmann constant, 1.3805 x 10 -23 J/K 
magnetic multipole of multipolarity l 
number of particles 

virtual photon number spectrum of a particular multipole, MeV -1 

neutron 

proton 

particular particle species emitted in heavy-ion collision 
alpha particle 

giant dipole resonance width, MeV 

maximum energy available for emission of secondary particles, MeV 
nuclear temperature 

particular electromagnetic multipole (e.g., El, MO) 
energy level density, MeV -1 

sum of cross sections in which at least one particle b is emitted 

photoneutron cross section, mb 

total photoneutron yield cross section, mb 

photonuclear total absorption cross section, mb 

total nucleus-nucleus EM absorption cross section, mb 

nucleus-nucleus EM reaction cross section for production of particular 
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nucleus-nucleus EM total absorption cross section summed over all EM 
multipoles, mb 
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^EM-abs nucleus-nucleus EM total absorption cross section for particular multipole, 
mb 

Cm giant dipole resonance cross-section peak value, mb 

<j* l (E,X) photonuclear reaction cross section for production of particular species X 
due to nuclear excitation from particular multipole nl, mb 
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Abstract 


The most important particle emission processes for electromagnetic 
excitations in nucleus-nucleus collisions are the ejection of single neu- 
trons and protons and also pairs of neutrons and protons. Methods are 
presented for calculating two-neutron emission cross sections in photo- 
nuclear reactions. The results are in a form suitable for application to 
nucleus-nucleus reactions. 


Introduction 

When cosmic rays in the form of heavy nuclei 
pass through spacecraft walls and astronauts’ bodies, 
they undergo an interaction with the atomic nuclei in 
the spacecraft or the astronauts. Of the four forces 
that we currently know about (strong, weak, electro- 
magnetic, and gravitational), the cosmic ray inter- 
action occurs via both the strong and electro- 
magnetic forces. The strong and electromagnetic 
(EM) cross sections are of comparable magnitude 
in some situations. Previous work (refs. 1 through 
5) has concentrated on studying this electromag- 
netic aspect of nucleus- nucleus collisions. This study 
has so far considered only single-nucleon emission 
processes. A preliminary study has been made of 
multiple-nucleon emission (ref. 5). However, the re- 
sults of this study are unable to be directly utilized 
in heavy-ion transport codes because experimental 
photonuclear cross sections are used as inputs into 
the calculation of nucleus-nucleus EM cross sections. 

The aim of the present work is to show how to cal- 
culate these photonuclear cross sections for multiple- 
nucleon emission. Given analytic expressions for 
these cross sections, it will then be possible to add 
multiple-nucleon emission due to the EM effect into 
heavy- ion transport codes. Calculating multiple- 
nucleon emission effects for the EM interaction is 
much more difficult than the calculation for single- 
nucleon emission. Thus the present work will con- 
sider only the most important multiple-nucleon emis- 
sion process— that of two-neutron emission. Other 
multiple-nucleon effects such as emission of two pro- 
tons, a neutron and a proton, or an alpha particle are 
more strongly suppressed than two-neutron emission 
primarily due to the Coulomb barrier. This is espe- 
cially true for neavy nuclei. In fact even single-proton 
emission is completely Coulomb-suppressed for heavy 
nuclei (ref. 1). 

Multiple-Nucleon Emission and Arbitrary 
Electromagnetic Multipoles 

Before proceeding to the study of two-neutron 
emission, we shall first place it in the broad con- 


text of arbitrary EM multipoles and multiple-nucleon 
emission. 

The total nucleus- nucleus EM absorption cross 
section <r EM is given by 

a EM = X^ EM ^ 

X 

where the sum is over all possible species emitted in 
the heavy- ion collison and ctem(^) the nucleus- 
nucleus EM reaction cross section for producing a 
particular species X. The cross section 0 em(^O 
was given as equation (1) in reference 1, but we 
now generalize it as (see also eq. (2.1) of ref. 6 and 
eq. (4.11) of ref. 7) 

= Y. L ,°* {E ' X) N "W dE ( 2 ) 

n l JE 0 \X) 

where tt can be either (E) or magnetic ( M ) and l is 
the order of multipolarity. Each term in the summa- 
tion of equation (2) represents nuclear excitation by 
a particular EM multipole. Each different EM multi- 
pole 7 rl causes a particular type of nuclear excitation. 
For example, the El photon field causes the nucleus 
to go into the giant electric dipole resonance (GDR) 
mode of oscillation parameterized as 

a ^ siE) = 1 + p2 - e£ dr ) /£ 2 r 2 ] (3) 

(see eq. (8) of ref. 1 and eq. (6) in this paper) 
which would subsequently decay into various chan- 
nels A, with a probability g Jp, sometimes called the 
branching ratio. However, an E2 multipole (ref. 8) 
would cause the nucleus to go into a giant electric 
quadrupole resonance (GQR), the cross section of 
which would have a different parameterization from 
equation (3), and where various decay probabilities 
g | 2 (e.g., neutron versus proton decay) may also be 
different from g jp. Thus the photonuclear reaction 
cross section for production of species X is some 
fraction of the total photonuclear absorption cross 
section 

a?(E,X) = g* x l (E) *iUe) 


(4) 


for a particular multipole. Equation (4) is a gener- 
alization of equation (7) of reference 1 and means 
physically that a particular EM multipole nl causes 
a collective nuclear vibration a^ s (E) which can then 
decay into various channels via g x (E). Note that 

£<$(£) = 1 ( 5 ) 

X 

By combining equations (2) and (4) and assuming 
an energy independent branching ratio (replacing 
g x (E) with g J?), the nucleus-nucleus EM reaction 
cross section becomes 

* £ $ [ <U E ) dE = £ 4'^' M . abs 
rt Je °W Tr 

( 6 ) 

where the nucleus- nucleus EM total absorption cross 
section is 

*EM-abs = / ^ <U E ) N*‘(E) dE (7) 
J hj 0 yX ) 

Equations (6) and (7) are generalizations of the equa- 
tions on page 7 of reference 1. (Note that in ref. 1 the 
absorption cross section was written as ^EM-abs(^) 
in equation (15). The dependence on X came out 
because the threshold E 0 (X) is the lower limit of in- 
tegration. We prefer here simply to write ^EM-abs 
because a true absorption cross section should not 
depend on X .) In previous equations, it is the photo- 
nuclear total absorption cross section cr^ g (£) which 
gets parameterized according to the particular nu- 
clear multipole excitation. For example, the nuclear 
El GDR excitation is parameterized in equation (3). 

If one assumes that the branching ratios are in- 
dependent of multipolarity (which may receive some 
justification from the Bohr independence hypothesis 
(ref. 9)), then, replacing g 1 ^ with g x , equation (6) 
becomes 

<r EM (£) « 9X £ °EM-abs = SX^EM-abs (8) 

nl 

where the nucleus-nucleus EM total absorption cross 
section summed over all EM multipoles is 

°EM-abs = £ °EM-abs ( 9 ) 

nl 

Equation (8) is interpreted physically as meaning 
that a nucleus-nucleus EM reaction has occurred 
which excites a nucleus into a superposition of multi- 
polarities ^EM-abs ( suc ^ as a Unear combination of 
nl 


GDR and GQR). This superposition can be likened 
to the compound nucleus concept of Blatt and Weis- 
skopf (ref. 9). This compound nucleus superposition 
then decays via multipole-independent branching ra- 
tios g x as in equation (8). 

Although equation (6) is probably more correct, 
it will be impossible to implement in practice due 
to the difficulty of calculating multipole-dependent 
branching ratios. Therefore in practical calculations, 
equation (8) will be used. This equation is entirely 
consistent with the concept of compound nucleus 
formation (ref. 9) and decay independent of the mode 
of formation (Bohr independence hypothesis). This 
is identical in spirit to the abrasion-ablation model 
(ref. 10). In fact, attempts (ref. 11) have been made 
to understand the abrasion-ablation model in terms 
of compound nucleus formation and decay using a 
T-matrix approach. Thus our basic equation (8) 
can be thought of as an electromagnetic abrasion- 
ablation model. 

Finally, we put our previous studies (refs. 1 
through 5) in the context of equations (1) and (8). 
First of all, equation (8) tells us that in order to cal- 
culate the nucleus- nucleus EM reaction cross section, 
we should sum over all possible EM multipolarities 
(£0, E 1, E 2, • • , MO, Ml, M2 , • * ■). The only multi- 
polarity that has been studied so far is El (refs. 1 
through 5) leading to the giant electric dipole res- 
onance. Clearly the effects of other multipolarities 
must be considered. Bertulani and Baur have con- 
cluded (ref. 6) that the electric quadrupole (E2) 
contribution to the total nucleus-nucleus EM cross 
section can be as much as 50 percent of the El con- 
tribution at 100 MeV/N and about 20 percent at 
1 GeV/N dropping to about 10 percent at higher en- 
ergies. Note that all the EM data (see ref. 1) are 
at high energy, and, thus, in comparing our theory 
(ref. 1) with experiment, the E2 contribution has 
not been large. Nevertheless we require the nucleus- 
nucleus (fragmentation and EM) theory to include 
all energies in the cosmic ray spectrum, and thus it 
is very important to consider other EM multipoles 
as well. Apart from mujtipoles other than El, equa- 
tion (1) tells us that we need also consider not only 
single- nucleon emission but multiple- nucleon emis- 
sion as well, such as emission of np , 2n, 2p, a, 
nnp , npp , and 3 n. However, references 1 through 
4 have considered only the £*1 multipole and only 
single-neutron emission. Thus, there still remains 
much territory to explore, namely multiple-nucleon 
emission with El excitation and then single- nucleon 
and multiple-nucleon emissions for all other multi- 
poles. The present paper is concerned only with 
two- neutron (2 n) emission from £1 excitation. 
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Theory of Two-Neutron Emission 

Cucinotta et al. (ref. 5) have in fact studied the 
problem of multiple-nucleon emission (2 n, 2 p, p3n, 
3 n) with El excitations only. However, the photo- 
nuclear reaction cross section < t v {E , X ) was obtained 
from experimental data, and even though useful con- 
clusions were drawn, the work of reference 5 was in- 
complete and cannot be used in a transport code, 
which requires analytical expressions for cr* l (E, X). 
The aim of the present work therefore is to begin 
the study of multiple- nucleon emission in nucleus- 
nucleus electric dipole (El) excitation reactions by 
using an analytic approach. The motivation for such 
an approach is to implement analytic cross-section 
expressions into cosmic ray transport codes. There 
is no cosmic ray transport code in existence which 
includes anything other than single-nucleon emission 
for EM reactions. Because a full theory for multiple- 
nucleon emission is much more complicated than for 
single- nucleon emission, the present work is limited 
to a study of two-neutron emission, which is the dom- 
inant multiple-neutron contribution. 

Two-Neutron Multiplicity 

Fuller et al. (ref. 12, p. 190; ref. 13, p. 4; ref. 14, 
pp. 143-145) have defined the total photoneutron 
yield cross section as 

<tH, Xn) = <7(7, n) + 2(7(7, 2n) + (7(7, np) + M7, 3n) + <7(7, n2p) + "' 

( 10 ) 

and the photoneutron cross section (i.e., the sum 
of cross sections in which at least one neutron is 
emitted) as 


the mulitplicity is 2 as expected. Often instead of 
using multiplicity, one works with the quantity 


<r( 7.2n) 

(7(7, n) + (7(7, 2 n) 


= M(E) - 1 


(14) 


Blatt and Weisskopf (ref. 9, ch. VIII, sec. 6B) have 
worked out the general theory of multiple-nucleon 
emission, which they call secondary nuclear reac- 
tions. They write the cross section for production 
of particle b (see ref. 9, eq. on p. 375), which we spe- 
cialize to photoproduction as 


<7(7, sb) = <7(7, b) + ^<7(7, 6c) (15) 

c 


This is made up of a primary, single-step cross sec- 
tion <7(7,6) in which particle b is emitted directly 
from decay of the compound nucleus and a secondary, 
multiple-step cross section <7(7,60) in which either 
particle 6 or c is emitted from the compound nucleus 
which decays to a lower excited state which again 
decays via emission of particle c or 6. Blatt and 
Weisskopf (ref. 9, pp. 376 379) then go on to cal- 
culate <7(7,6) and <7(7,60) in terms of <7(7, s6); how- 
ever, one generally needs to work the integrals out 
numerically. We shall eventually do this when work- 
ing out the general expression for multiple-nucleon 
emission. However, for the case of two-neutron emis- 
sion, certain simplifying assumptions concerning the 
integrals can be made so that analytic expressions 
can be obtained. These are (ref. 9, eq. (6.14), p. 377) 


< 7 ( 7 , n) «<7(7,sn)(l + ^)exp(-|^) (16) 


cr(7, sn) = <7(7, n) 4- <7(7, np) 4* <?( 7> 2rt) + a{j, na) 4- <7(7, 3n) + • ** 
The neutron multiplicity is defined (ref. 15) as 


M(E) ^ 


( 12 ) 


<7(7, sn) 

which for n and 2 n emission only becomes (ref. 16) 


and (ref. 9, eq. (6.18), p. 379) 

<7(7, 2n) ~ <7(7, sn) 1 - ^1 + exp ^ 


where © is the nuclear temperature and 
fsec = E'y Egipj, 2n) 


(17) 

(18) 


M{E) « 


<7(7, n) + 2<7(7, 2n) 
<7(7, n) -|- <7(7, 2n) 


(13) 


This expression for the multiplicity has been implic- 
itly assumed also in Thies and Spicer (ref. 17, eq. (5)). 
Note the reason for writing the neutron multiplicity 
this way: for <7(7, 2n) equal to zero (meaning that 
only single neutrons (n) are emitted), the multiplic- 
ity is 1 as expected, and for <7(7,71) equal to zero 
(meaning only two- neutron pairs (2n) are emitted), 


Equations (16), (17), and (18) are also discussed in 
references 15 through 18. Equation (17) is derived 
on the assumption that the photoneutrons are pro- 
duced via statistical decay of the compound nucleus 
(ref. 9). However this is not always necessarily true. 
It can sometimes happen that the incident photon, 
rather than exciting a compound nucleus which sub- 
sequently decays via neutron emission, knocks a neu- 
tron out directly. (Note that this direct emission vi- 
olates the abrasion-ablation concept.) We therefore 
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introduce a fraction of direct emission fa, which typ- 
ically has values (ref. 15 and private communication 
from G. O’Keefe and R. Rassool of the University of 
Melbourne) in the range 0. 1-0.2 although it can be as 
large as 0.4 (ref. 17). This direct fraction is incorpo- 
rated into our theory (refs. 15 and 17) by generalizing 
equation (17) as 

ff (7.2«)*a(7.»0(l-,fc) [l - (l + exp )] (19) 

This equation is our main result for 2 n emission. The 
quantities fa, e sec , and 0 are easily calculated so 
that the only input required is <7(7 ,sn). Blatt and 
Weisskopf (ref. 9, p. 379) claim that (7(7, sn) can 
usually be closely approximated by the compound 
nucleus formation cross section, that is, the total 
absorption cross section. (This is done in ref. 18.) 
However, this is only true when compound nucleus 
decay proceeds predominantly by neutron emission 
which is the case for heavy nuclei. For light nuclei, 
where proton emission is just as important (ref. 1), 
this approximation is not valid. Thus we now neglect 
all contributions to a(j , sn) other than <7(7,71) and 
<7(7, 2n) so that 


^(7>2n) 

<t( 7, n) + <7(7, 2n) 


w (1 - fd) 



t (20) 
which is the result written in references 16 and 18 for 
fd = 0. Defining 


Xfrac = (1 - f d ) [l - (l + *S) exp (-ip Y 

( 21 ) 

we rearrange equation (20) to obtain 

, - . ATrac , , 

<7(7 ' 2 " )= l-Xfrac"* 7 ’") (22> 

Although equation (19) is our fundamental equation 
for 2n emission, we shall in fact use equation (22) to 
test the general theory. This is because several au- 
thors have reported both <7(7, n) and <7(7,2 n). We 
shall use experimental input for <7(7, n) and then cal- 
culate <7(7, 2n) by equation (22) and see how well our 
calculations agree with experimental measurements. 
If we obtain good results, we can be confident that 
the theory outlined can be used to calculate 2n emis- 
sion in EM nucleus- nucleus collisions. 

Finally, we note that equations (21) and (22) can 
be used for an approximate calculation of the neutron 


multiplicity given in equation (13) as 


mm = 1 + a - s d ) [1 - (1 + e -f) exp (-|p)' 

, . . (23) 

which is the result reported in references 15 and 17, 
(Note however the error in eq. (3.12) of reference 15, 
which is that <7(7, np) should not appear in the 
numerator.) 

Nuclear Temperature 

In the previous sections, the nuclear temperature 
0 plays a central role. For a Maxwell-Boltzmann 
(classical) gas (ref. 19, p. 117), the relation between 
energy and temperature is linear 


E = ^N 0 kQ (24) 

whereas for a Fermi gas, the relation is quadratic 

E — (Constant)& 2 0 2 (25) 

The theory developed by Blatt and Weisskopf (ref. 9, 
p. 372) uses a Fermi gas nuclear temperature given 
by 

_ fE \ 1/2 

® = ( 1 ) < 26 > 

where the constant A depends on the total nucleon 
number A. We shall not use the Blatt and Weisskopf 
theory for nuclear temperature, but rather use a 
more sophisticated version developed by Bohr and 
Mottelson where (ref. 20, eq. (2-50), p. 154) 


© — 1 


1 dp 

pdE 


(27) 


with the Fermi gas energy level density (ref. 20, 
eq. (2-47), p. 153) 


P(N, Z, E) 


6 1 /* go 
12 ( So E)5/4 6XP 



(JVwZ) 


(28) 

being of similar form to that of Blatt and Weisskopf 
(ref. 9, eq. (2-50), p. 154). The symbol E is the ex- 
citation energy, N is the neutron number, Z is the 
proton number, and g 0 is the one-particle level den- 
sity at the Fermi energy, representing the sum of neu- 
tron and proton level densities (ref. 20, p. 153). Thus 
the nuclear temperature becomes (ref. 20, eq. (2-50), 
p. 154) 



- 7r (29) 
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which is the same form as used in reference 15. It 
is very important to realize that the energy E in 
equations (28) and (29) is, in fact 

E = E 1 - E 0 (rf,n) (30) 

given by the single-neutron emission threshold 
E 0 { 7, n)- This is emphasized in reference 18 (p. 304) 
and by G. O’Keefe and R. Rassool (private communi- 
cation). The one-particle level density g 0 is related to 
the level density parameter a via (ref. 20, eq. (2-123), 
p. 187) 



which is plotted in figure 2-12 of reference 20 (p. 187). 
It is seen that a very much depends on shell struc- 
ture, and therefore whenever possible, a should be 
obtained directly from figure 2-12 of reference 20. 
However, a rough approximation is 


Table 1. 12 C Photoneutron Cross Sections 
[Data are from ref. 21] 


MeV 

7.«)» 

mb 

< 7 ( 7 , 2 n), mb 

Experimental 

Theoretical 3 

35 

2 

0.004 

0.04 

40 

.8 

.012 

.05 

50 

.2 

.015 

.02 

60 

.1 

.014 

.01 

80 

.1 

.095 

.01 


a fd — 0.9; a = 2 MeV -1 


Table 2. 18 O Photoneutron Cross Sections 
[Data are from ref. 22] 


am- MeV -1 (32) 

8 

but note that this fails badly for A between 190 
and 210 (ref. 20). When the nuclear temperature 
is rewritten as 


\eJ 4 E 


(33) 


with the first term corresponding to the form used 
by Blatt and Weisskopf (ref. 9), we approximate a to 
obtain 


1 - fAY /2 _ A 


6 U EJ 


4E 


(34) 


E'y, 

MeV 

< 7 ( 7 , n), 
mb 

<y( 7 , 2 n), mb 

Experimental 

Theoretical 3 

14 

4 

1 

0.5 

18 

2.4 

3 

1 

22 

4 

3 

3 

26 

7 

3 

6 

30 

5 

3 

5 

34 

2 

3 

2 

38 

2 

2 

2 

42 

1 

3 

1 


°/j = 0.5; o = 3 MeV -1 


Test of Two-Neutron Multiplicity Theory 

A computer code was written to implement the 
theory described. Necessary inputs are the energy 
level density parameter a which was taken directly 
from figure 2-12 (ref. 20, p. 187) rather than from 
equation (32). The two-neutron cross section was 
then calculated from the single-neutron cross section 
with a best-fit value of f d , the fraction of direct 
emission. The results are listed in tables 1 through 
8. (The data in tables 1 through 8 are taken from 
references 16 and 21 through 27.) As can be seen, the 
agreement between theory and experiment is good, 
except for a few notable exceptions. It is disturbing, 
however, that the direct emission fractions are so 
high. 


Table 3. 40 Ar Photoneutron Cross Sections 
[Data are from ref. 23] 


Ey, 

MeV 

<7(7, *0. 

mb 

< 7 ( 7 , 2 n), mb 

Experimental 

Theoretical 3 

18 

16 

5 

6 

20 

7 

14 

12 

21 

5 

13 

14 

22 

3 

11 

12 

23 

2 

10 

11 


°/<f = 0.05;o = 6MeV“ 1 
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Table 4. 59 Co Photoneutron Cross Sections 
[Data are from ref. 24] 


E*y, 

MeV 

<7(7, n), 
mb 

<7(7, 2n), mb 

Experimental 

Theoretical 0 


47 

9 

4 


21 

17 

7 


17 

12 

9 


14 

11 

8 


19 

9 

12 


20 

5 

13 


9 

6 

6 


1 

6 

1 


4 

3 

3 


% = 0.6; a = 8 MeV _1 


Table 5. ^Cu Photoneutron Cross Sections 
[Data are from ref. 25] 


E 7, 
MeV 

cr(7,n), 

mb 

<7(7, 2n), mb 

Experimental 

Theoretical 0 

20 

44 

2 

0.6 

21 

26 

4 

5 

22 

20 

8 

11 

23 

13 

12 

12 

24 

11 

13 

13 

25 

10 

13 

15 


% = 0.3;a = 8 MeV- 1 


Table 6. ® 9 Y Photoneutron Cross Sections 
[Data are from ref. 26] 


E'y, 

MeV 

<7(7, n), 
mb 

<7(7,2 n), mb 

Experimental 

Theoretical 0 

22 

35 

7 

6 

24 

10 

10 

6 

26 

2 

16 

2 

28 

.5 

16 

.5 


a fd — 0*5; a = 12 MeV* 1 


Table 7. 197 Au Photoneutron Cross Sections 

Data are from ref. 16 but have been multiplied by 0.93 
following the suggestion of ref. 27 


MeV 

<7(7, n), 
mb 

<7(7, 2n), mb 

Experimental 

Theoretical 0 

15 

391 

9 

20 

16 

214 

73 

121 

17 

93 

98 

98 

18 

56 

93 

74 

19 

37 

65 

53 

20 

27 

58 

40 

21 

17 

42 

25 

22 

11 

47 

16 

23 

7 

40 

11 

24 

2 

37 

3 


a fd = 0.4; a = 24 MeV -1 


Table 8. 208 Pb Photoneutron Cross Sections 
[Data are from ref. 16] 


Ei, 

MeV 

<7(7, «), 

mb 

<7(7, 2rc), mb 

Experimental 

Theoretical 0 

16 

180 

80 

104 

18 

30 

80 

51 

20 

20 

50 

55 


a f d = 0.2; a = 8 MeV -1 


Discussion 

We have demonstrated that the theory of two- 
neutron emission described herein is able to predict 
most two-neutron multiplicities with a high degree 
of accuracy. One therefore wishes to implement this 
theory into the calculation of nucleus-nucleus cross 
sections and eventually into a transport code. A pa- 
rameterization of the energy level density parameter 
a can be very easily worked out (ref. 20 ). Then the 
only unknowns are and <7(7, n). 

The values of were obtained from an overall 
best fit to <7(7, 2 n) data. (Ideally, a single <7(7, 2 n) 
value at a single energy determines f d .) Clearly some 
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further work remains to be done here. One can either 
work out a parameterization by determining fa from 
best fits to a whole range of data, or better still, one 
should determine fa theoretically. Our initial effort 
will undoubtedly be a parameterization so that the 
theory can be used within transport codes. 

The most serious unknown is <7(7, n) which is 
used as input for the (7(7, 2n) calculation at each en- 
ergy. One might think that the liquid drop model 
calculations for (7(7,71), described extensively in ref- 
erence 1, would be adequate for this purpose, es- 
pecially as these same model calculations gave such 
good results for nucleus-nucleus cross sections. Un- 
fortunately, this is not the case. Even though the 
liquid drop model calculations of (7(7, n) are quite 
accurate at most photon energies, they are unfortu- 
nately not very accurate for high photon energies. 
The reason that this did not matter for the calcu- 
lations of reference 1 was that the virtual photon 
spectrum N(E) predominated at low photon ener- 
gies, and therefore energy integrated nucleus-nucleus 
cross sections were quite accurate. However, obvi- 
ously (7, 2n) processes only occur in the high energy 
(7,n) region. The calculations of <7(7, 2 n) require 
<7(7, n) as input, and thus we now need very accu- 
rate calculations for (7(7, n) at high energy. (Calcu- 
lations described in the present paper used experi- 
mental numbers.) 

Thus, the most pressing need for inclusion of 
the effects of multiple-nucleon emission in heavy-ion 
transport codes is the accurate calculation of (7(7, n) 
at high energy. The first such step will require the 
study of deformation splitting (refs. 14 , 17 , 18 , 28 , 


and 29 ) and isospin splitting of the giant dipole 
resonance (refs. 28 and 29 ). 

Finally, even though only two-neutron emission 
has been considered herein, it is now clear how to 
proceed with the general multiple-nucleon emission 
problem. One can simply numerically integrate the 
expressions of Blatt and Weisskopf (ref. 9 ) to obtain 
any multiple-nucleon final state. (These expressions 
simplified in the two-neutron case considered herein, 
so that numerical integration was not necessary.) 
This method has previously been applied to 32 5(7, d) 
by Norbury (ref. 30 ). 

As an alternative to the whole approach taken 
in the present work (and even for the single- nucleon 
work of ref. 1 ), it may be feasible to instead do calcu- 
lations using some of the modern nuclear evaporation 
codes (e.g., EVAP -4 and EVA- 3 ). The advantage of 
these codes is that they fully incorporate both statis- 
tical and direct emissions. The excitation energy in 
these codes is simply the photon energy considered 
above. 

Concluding Remarks 

Electromagnetic excitations in nucleus- nucleus 
collisions can often involve the ejection of pairs of nu- 
cleons. Methods have been presented for calculating 
two-neutron emission cross sections in photonuclear 
reactions. These cross sections are now in a form 
suitable for use in nucleus- nucleus reaction theory. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November 20, 1989 
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